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Abstract. Measurements of inclusive hadron suppression and di-hadron azimuthal 
correlations in ultra-relativistic nuclear collisions at RHIC have provided impor- 
tant insights into jet quenching in hot QCD matter, but are limited in their sen- 
sitivity due to well-known biases. Full jet reconstruction in heavy-ion collisions 
would conceptually provide a direct measurement of the energy of the scattered 
parton before energy loss, alleviating such biases and allowing a measurement of 
^ ' the energy loss probability distribution in a model-independent way from hard 

probes. In these proceedings we utilize recent progress in the reconstruction of 
jets in the heavy ion environment and present the first measurement of the frag- 
mentation function from fully reconstructed jets in heavy ion collisions. The frag- 
^ , mentation function measured in central Au+Au collisions at 200 GeV will 

be presented and discussed with respect to p+p reference measurements. 
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1 Introduction 



Measurements of inclusive hadron suppression and di-hadron azimuthal correlations in ultra- 
relativistic nuclear collisions at RHIC have provided important insights into the properties of 
hot QCD matter [1 2j. In particular, the high transverse momentum (high pt) suppression [3] 
and low pt enhancement [4J of the correlated yield of hadrons recoiling from a high pt particle 
\ (azimuthal pair separation Acj) ~ tt) suggest a dramatic softening of jet fragmentation in dense 

. matter, arising from strong partonic energy loss. However, these measurements are limited in 

■ their sensitivity due to well-know geometric biases (see for example [5j). Full jet reconstruc- 

tion in heavy-ion collisions would conceptually provide a direct measurement of the energy of 
the scattered parton before energy loss, alleviating such biases and allowing to reconstruct the 
partonic kinematics in an unbiased way, independent of the fragmentation details (quenched or 
unquenched). This would result in a direct measurement of the energy loss probability distri- 
bution necessary in a model-independent way. 

A sensitive new observable accessible via full jet reconstruction is the fragmentation function 
and it's modification in the presence of jet-quenching. The fragmentation function expressed 
in the logarithmic variable ^ = \n{l/x) = \n{Ejet/p^) with x = p^/Ejet^ commonly known 
as the "hump-backed" plateau, is well-described by QCD-based MLLA (modified leading-log 
approximation) calculation in e+e~ collisions [6j. Due to current experimental constraints (see 
sec. 12.2. ip we will express the ^ fragmentation function using 

e = ln(pI-,/p^'^'-°"). (1) 



^ e-mail: joerii.putschke@yale.edu 
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Fig. 1. (color online) Single event display of the summed charged and neutral pt per grid cell in the 
?7,(/) plane from 0-20% central HT triggered Au+Au collisions. Left panel: reconstructed di-jet with 
Pir/et- 21 GeV. Right panel: single jet with pl%t^ 47 GeV. 



Jet quenching in the MLLA formalism has been introduced via modification of the radiation 
vertex coupling strength [6j. In these calculations a significant enhancement in the fragmenta- 
tion function at high ^ (low hadron pt) as well as a suppression at low ^ (high hadron pt) with 
respect to p+p (or e+e~) reference measurements is expected. 

In these proceedings we will report the first fragmentation function measurement from full 
jet reconstruction in central Au+Au collisions at ^/s^= 200 GeV. The heavy- ion fragmen- 
tation function will be compared and discussed with respect to p+p reference measurements. 
In addition we will discuss the effect of the underlying high multiplicity heavy-ion background 
on the jet energy determination as well as appropriate background subtraction methods in 
heavy-ion collisions to measure the jet energy and the fragmentation function. See [7J for the 
accompanying discussion about inclusive jet spectra and systematic studies by utilizing differ- 
ent jet reconstruction algorithms in heavy-ion collisions, also presented for the first time during 
the Hard Probes 2008 meeting. 



2 Experimental setup and jet reconstruction performance 
2.1 Experimental setup and event selection 

The full jet reconstruction measurements presented in these proceedings were performed by 
combining the charged particles from the STAR TPC and the neutral (tower, 0.05 x 0.05 in 
T] X (j)) energy from the EMCal. Both detectors, TPC and EMCal, have full azimuthal coverage 
(0 < ^ < 27r) and an acceptance at mid-rapidity ranging from — 1 < < 1 (see Fig. [T]). TPC 
primary- vertex tracks within < 1 were selected for this analysis using standard quality cuts 
for high-p^ measurements [2] , which eliminate fake tracks and ensure sufficient momentum reso- 
lution. For the neutral tower energy corrections were applied to account for energy deposited by 
hadrons as well as to avoid double-counting of the electron energy. For this correction electrons 
were identified using d.p/E criterium and only the pt of the TPC electron track was used further 
on. In addition these identified electrons were not used in the charged particle fragmentation 
function measurements, because the majority of these are from ttq conversions in the STAR 
detector material. 



The following event sets were used for this analysis, based on the online trigger conditions: 
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Fig. 2. (color online) Left panel: Background fluctuations in 0-20% central MB Au+Au collisions 
as function of the cone radius R for different pt^^ (more details see text). Right panel: Reconstructed 
Pythia jet spectrum with and without embedding in 0-20% central MB Au+Au collisions. 

— (MB Au+Au) minimum-bias trigger in Au+Au collisions: year 7 RHIC run coincidence 
condition between the two ZDCs in STAR. 

— (HT Au+Au) HT-trigger in Au+Au collisions: year 7 RHIC run with the trigger condi- 
tion that the sum of the two highest towers in a 3x3 patch (with Et > 5.67 GeV for the 
highest tower) satisfy Et > 7.5 GeV (in addition to the minimum-bias trigger condition); 
sampling an integrated luminosity of ~ 500 which corresponds to a p+p equivalent of 
- 19.6 pb-^. 

— (HT p+p) HT-trigger in p+p collisions: year 6 RHIC run with the trigger condition for 
a single EMCal tower Ef > 5.4 GeV (in addition to the minimum-bias trigger condition), 
corresponding to an integrated luminosity of ~ 11 pb~^. 

2.2 Jet reconstruction algorithms and performance 

2.2.1 Jet reconstruction algorithms 

If not otherwise mentioned jet-finding was performed using a seeded cone algorithm which 
combines all particles in the r] x (j) plane within a cone radius of Rc = \/A^^+Arf around a 
high energy seed. In the presence of heavy-ion background splitting/merging and the iteration 
procedure to define stable cones can lead to ambiguities. As a first step we utilized a cone 
algorithm without splitting/merging and iteration: leading order (high seed) cone (LOConefl 
The relevant parameters for this algorithm are and Rc- If not otherwise mentioned a 

pseed_ ^ Q GeV and Rc= 0.4 was used for all measurements. In addition a conceptually dif- 
ferent algorithm, based on sequential recombination, was used. We choose the kt algorithm, 
part of the Fast Jet package (for more details see [8] and references therein). For the LOCone 
a grid, based on the EMCal tower granularity (0.05 x 0.05 in x 0), was filled with charged 
particles and the neutral tower energy. Jet finding was performed on that grid with the summed 
Pt per grid celQ. See Fig. [T] for two 0-20% central HT triggered Au+Au events containing jets 
reconstructed using the LOCone. 

In these proceedings only the jet with the highest reconstructed jet energy per event (in 
simulations and data) was used for jet spectra and fragmentation function measurements. 

^ http : //rhig. physics. yale.edu/ ~ putschke/Ahijf / 

^ The Pt from charged particles was used (no PID) and Et = pt for the neutral tower energy, assuming 
that a photon deposited the tower energy. 
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Fig. 3. (color online) Left panel: ^ fragmentation function for 30 GeV Pythia jets embedded in 0-20% 
central MB Au+Au collisions. The red dashed line represents the background contribution (for more 
details see text). Right panel: ^ fragmentation function for 30 GeV Pythia jets embedded in 0-20% 
central MB Au+Au collisions after background subtraction compared to the Pythia fragmentation 
function (for reference the non-background subtracted fragmentation function is show as the dashed 
histogram) . 



2.2.2 Heavy-ion background and effect on jet spectrum 

The main complication in performing full jet reconstruction in heavy- ion collision is the presence 
of the large underlying background. The common approach is to reduce the size of the cone 
radius Rc to 0.3 — 0.4 instead of 0.7 as normally used in p+p collisions. This reduction will 
still recover ~ 80% of the jet energy within Rc ~ 0.3 for 50 GeV jets as observed in p+p 
collisions at the Tevatron [9j. Jet quenching in heavy- ion collision is expected to broaden the 
jet fragmentation and would reduce the fraction of the measured jet energy at a given cone size. 
This effect needs further experimental studies and theoretical exploration. To further reduce 
the heavy- ion background a minimum pt cut (p^^^) of up to 2 GeV on charged particles and 
neutral towers is applied jlO|llj . The remaining mean background (< pf^^' >) is estimated 
in the out-of-cone area, excluding the area occupied by the two highest reconstructed jets per 
evenl0, and corrected for by scaling to the proper jet cone area. 

To study in more detail how the underlying heavy-ion background correction depends on 
Rc and p^'^^ we measured the region-to-region background fluctuations in 0-20% central MB 
Au+Au events. The background fluctuation measurement is characterized by the gaussian 
width: a{Apt{Rc,pr*)) = ^^'''^^^(i^cPr^))- < Pt^^' >• In Fig.[2]left panel the background 
fluctuation for different p^'^^ is shown as function of Rc- One can see that applying a p^^* of 
2 GeV significantly reduces the background fluctuation from ~ 6 GeV for Rc = 0.4 with p^^^ 
= 0.1 GeV to below 2 GeV for p^^^^ = 2 GeV . 

The effect of background fluctuations on the jet energy resolution and on the reconstructed 
jet spectrum was studied by embedding a Pythia generated jet-spectrum [12 {Pythia) into 
0-20% central MB Au+Au collisions {Pythia+AuAu 0-20%). For the jet signal the charged 
and electromagnetic particles were taken directly form PythiaEI. In Fig. [2] right panel the re- 
constructed Pythia jet spectrum with and without embedding in 0-20% central MB Au+Au 
collisions is shown for a p^^^= 0.1 and 2 GeV. The jet energy resolution in heavy- ion collisions 
is mainly dominated by the background fluctuations (see [7J for more details). It is apparent 
that for a pl^^= 0.1 GeV the steeply falling jet spectrum folded with the jet energy resolu- 
tion — determined by the background fluctuations — leads to a significant broadening of the 

^ Strong elliptic flow V2 can lead to systematic effects in determining < pf^^' > depending on the 
relative location of the jet-axis with respect to the reaction plane. With sufficient acceptance and 
focusing on central Au+Au collisions, where V2 is small, the effect is expected to be small. 

^ The detector effects are expected to be small especially compared to the uncertainties caused by 
the heavy-ion background. 
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Fig. 4. (color online) Left panel: ^ fragmentation function in 0-10% central HT Au+Au collisions after 
background subtraction applying a Pt^*> 2 GeV for pl^jet> 30 GeV using the LOCone and Fast Jet 
kt jet-finding algorithm. The arrow indicates the ^ value for a hadron with pt^ 10 GeV. Right panel: 
^ fragmentation function in 0-20% central HT Au+Au collisions compared to HT p+p collisions for 
Pt%t{pp) > 30 GeV using the LOCone with Rc= 0.4 and a pr*> 2 GeV (for reference the non- 
background subtracted fragmentation function is show as the dashed histogram). The grey shaded 
area in both panels indicates the region in ^ (corresponding to p^^^^^^ < 2 GeV) where the Au+Au 
measurement is expected to be dominated by the underlying background (see discussion in sec. 12.2.3]) . 



reconstructed jet spectrum in Au+Au collisions compared to the input Pythia spectrum. As 
a consequence, for a pf^^= 0.1 GeV, jets with a reconstructed jet energy of pl^f^t^ 30 GeV in 
the presence of background will be dominated by the feed-down of lower Pt^fet J^^^- -"-^ contrast 
the jet spectra for a pf^^= 2 GeV with and without embedding in 0-20% central MB Au+Au 
collisions are in very good agreement, therefore allowing a good determination of p^j^e^ 
event-by-event basis in Pythia+AuAu 0-20% as compared to Pythia. 



2.2.3 Fragmentation function measurement performance 



To measure the ^ fragmentation function the jet energy pl^j^^t t>e determined on an 

event-by-event basis. To achieve this, one has to apply a pf^^ of 2 GeV to suppress the effect of 
background fluctuations, as discussed in sec. I2.2.2[ In the following the fragmentation function 
measurements were performed using the LOCone with Rc= 0.4 and the kt algorithm with 
Rparam = 0.4 and a pI^^= 2 GeV to determine the jet energy Pt^f^f contrast, all charged 
particles in a cone of R^^ = 0.7 (around the jet axis) without applying an additional pt 
were included in the fragmentation function measurements. In Fig. [3] left panel the charged ^ 
fragmentation function of a 30 GeV Pythia jet embedded in 0-20% central MB Au+Au collisions 
is shown. For ^ values larger then two the fragmentation function is dominated by the soft 
underlying heavy-ion background. The red dashed curve represents the ^ background estimate 
determined in the out-of-cone area (see sec. I2.2.2p on an event-by-event basis after scaling to 
pl^fet fragmentation function area. Fig. [3] right panel shows the comparison of the 

the background corrected ^ fragmentation function in Pythia+AuAu 0-20% to Pythia. The 
used background subtraction method and the jet energy resolution causes deviations of less 
then 10 — 20% for ^ < 2.5 assuming Pythia fragmentation. The deviations in the ^ shape for 
Pythia+AuAu 0-20% can be largely explained by jet energy resolution effects. For ^ > 2.5 — 3 
the fragmentation function measurement is dominated by the underlying heavy-ion background. 
Without further studies we will exclude these higher ^ regions in the following discussion. 
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Fig. 5. (color online) Ratio of the ^ fragmentation function measured in 0-20% central HT Au+Au 
collisions to HT p+p collisions for Pt^fetiPP) > 30 GeV with pt^^> 1 and 2 GeV (more details see text). 
The grey shaded area indicates the region in ^ (corresponding to p^^^^^^ < 2 GeV) where the Au+Au 
measurement is expected to be dominated by the underlying background (see discussion in sec. 12.2.3]) . 

3 Fragmentation function measurements in Au+Au (and p+p) collisions 

Using the algorithms and background correction procedure as discussed in sec. 12. 2| we per- 
formed the first fragmentation function measurement from full jet reconstruction in heavy-ion 
collisions using central HT Au+Au collisions at y/s^= 200 GeV. Fig. [H left panel shows the ^ 
fragmentation function for pl^f^t^ GeV in 0-10% central HT Au+Au collisions after back- 
ground subtraction using the LOCone and Fast Jet kt jet-finding algorithm. A pl^^> 2 GeV 
was applied and for the LOCone a Rc= 0.4 and for the FastJet kt a Rparam = 0.4 was used 
to suppress the effect of background fluctuations (see discussion in sec. 12.2.3]) . The comparison 
of the measured fragmentation function in 0-10% central HT Au+Au collisions for ^ < 2.5 by 
utilizing two conceptually different jet finding algorithm, LOCone and FastJet kt, given the 
applied parameters, is a first step to address systematic uncertainties in full jet reconstruction 
in heavy-ion collisions. Further systematic studies are needed, but the agreement of LOCone 
and FastJet kt gives us confidence in our jet reconstruction performed on central HT Au+Au 
collisions. 

In Fig. m right panel the ^ fragmentation function in 0-20% central HT Au+Au collisions 
compared to HT p+p collisions for Pt^f^tiPP) ^ GeV using the LOCone with Rc= 0.4 and a 
pcut^ 2 GeV is shown. In order to select a corresponding p+p jet population in 0-20% central 
HT Au+Au collisions Pythia simulations (see sec. 12. 2. 2j) were used to determine the appropriate 
pl^f^^{AuAu) region: pl^f^^{AuAu) > 31 GeV for p^''*> 2 GeV and pl^f^^{AuAu) > 35 GeV for 
pf^^> 1 GeV. 

To quantify in more detail the modifications in the fragmentation function due to quenching 
effects in Au+Au with respect to p+p reference measurements, the ^ fragmentation function 
ratio of Au+Au to p+p for pl^f^tiPP) ^ GeV is shown in Fig. [5l No apparent modification 
in the Au+Au ^ fragmentation function is seen in the ^ rang ^ 0.5 < ^ < 2.75 for p^^^^> 2 GeV 
(and 1 GeV). 



^ The lower ^ cut off is currently motivated by the momentum resolution of charged tracks above 
pt> 10 GeV. Further studies are needed to estimate and correct for the momentum resolution effects 
at highest pt. 
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Fig. 6. (color online) Inclusive jet spectrum for 0-10% MB and HT Au+Au collisions compared to Nhin 
scaled p+p reference measurements [13 . The MB Au+Au inclusive jet spectrum is correct for missing 
energy caused by applying a p^^* and Rc cut as well as for efficiency using Pythia simulations (for more 
details see text and [7 ). The HT Au+Au inclusive jet spectrum is not corrected. The orange shaded 
area indicates the p[jet region in which the ^ fragmentation function measurements were performed 
(see sec. [3}. 



One should note that the presented ^ distributions are not corrected for the missing jet 
energy imposed by applying a pf^^ (and Rc) cut. This fraction could be estimated assuming 
Pythia fragmentation (see sec. I2.2.2|) . but is expected to be different due to quenching effects 
in Au+Au collisions (more discussion in sec. [4]). 



4 Interpretation and discussion 

Assuming that full jet reconstruction allows us to select an unbiased jet sample {R'aX— ^) 
with the measured inclusive hadron suppression in heavy- ion collisions {Raa < 1) one would 
expect to the see modifications in the fragmentation function measurement due to jet quenching 
in heavy- ion collisions with respect to p+p reference measurements. In this light, our finding 
that there seems to be no modification in the fragmentation function in Au+Au collisions (see 
sec. [3]) is very surprising. To investigate further whether an unbiased sample of jets is selected 
in the fragmentation function analysis, we compare the jet spectra for 0-10% MB (corrected 
for p^^^, Rc and efficiency using Pythia simulations) and HT Au+Au collisions (no correction 
applied) to the TV^^^ scaled p+p reference measurements [13^ as shown in Fig. [6l Using a p^^^of 
2 GeV in the MB Au+Au jet spectrum measurement suppresses the jet yields with respect to 
the Ni)in scaled p+p jet spectrum. Applying a p^^^ in heavy- ion collisions leads to R'^a^ 
indicating a bias in the selected jet population in Au+Au collisions. However, measurements 
of the inclusive jet spectrum without applying an additional p^'^* indicate that we might be 
able to recover an unbiased jet sample {Raa^ 1) (see [7 ). This result suggests that in the 
presence of jet quenching in heavy-ion collision the correction for the p^'^^ and reduced Rc 
based on Pythia is not sufficient and does seem to underestimate the reconstructed jet energy 
in heavy- ion collisions. To illustrate this effect on the fragmentation function measurements 
see the sketch in Fig. [71 In the presence of large quenching correcting the jet energy based on 
unmodified jet fragmentation from Pythia will underestimate the reconstructed jet energy in 
central Au+Au collisions. The true jet energy in Au+Au is higher than in p+p collisions and 
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Fig. 7. (color online) Sketch to illustrate the effect of jet-energy bias due to energy loss and the 
applied pt^* on the measured ^ fragmentation function. 



will therefore lead to a relative shift of the ^ fragmentation function in Au+Au with respect 
to p+p, as illustrated in Fig. [3 right panel. This relative shift would then lead to the expected 
modification pattern of the fragmentation function in Au+Au due to jet quenching (see sec.[T]). 

An alternative scenario that would describe the data is the "black-and-white picture" , where 
some partons lose so much energy that they are "lost" (i^^^< 1), while the remaining fraction 
of jets suffers little or no modification, so that the fragmentation function of the "surviving" 
jets is unmodified. To distinguish between these two scenarios one has to establish if R'^X using 
our jet reconstruction framework is unity (or below). 

In addition the HT trigger requirement in p+p and Au+Au will lead to a trigger bias in the 
jet spectra and fragmentation functions. Recent studies in p+p collisions at ^/s^= 200 GeV 
[14] indicate that for a HT trigger with a tower energy threshold greater than 5.4 GeV biases 
in the fragmentation function are still present at pl^f^t {pp) ~ 30 GeV. To study the HT trigger 
bias in the presence of jet quenching effect, by even applying a higher trigger threshold of 7.5 
GeV, further statistics in the MB Au+Au data sample is needed (and available on tape) to 
assess the bias at higher pl^^et- "^^^^ trigger bias could lead — in the pl^fet I'ange of > 30 GeV 
used for the fragmentation function measurements — to a bias in the selected jet population 
favoring jets with no or only small energy loss. 



5 Summary and Outlook 

In these proceedings the first fragmentation function measurement from full jet reconstruction 
in heavy- ion collisions was presented (see sec. [3]). It was shown that using a smaller cone ra- 
dius of Rc= 0.4 and applying a p^^^of 2 GeV reduces significantly the effect of the underlying 
background fluctuations in our jet measurements, allowing — based on Pythia simulations (see 
sec. I2.2.2P — a reasonably precise determination of plj^^t an event-by-event basis necessary 
for the measurement of the fragmentation function. Utilizing the current jet finding algorithms 
and background subtraction procedures deviations of less then 10 — 20% for ^ < 2.5 assuming 
Pythia fragmentation are observed in the ^ fragmentation function measurement (see sec. 12.2.3]) . 
However, applying a p^'^^to suppress the heavy-ion background leads to a bias in the selected 
jet-population which would result in an underestimation of the true jet energy pl^-^^^{AuAu) in 
the presence of large jet quenching effects (see discussion in sec. |4]). This could explain why 
in our measurement we do not see an apparent modification in the fragmentation function 
in Au+Au collisions. To correct properly for this effect further experimental and theoretical 
developments are needed. First calculations addressing the effect of experimental cuts on the 
reconstructed jet energy in the presence of quenching can be found in [15j. 
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Fig. 8. (color online) Di-jet Acj) distribution for pl%t> 20 GeV in 0-20% central HT Au+Au and HT 
p+p collisions using the LOCone algorithm with Rc= 0.4, pr*= 2 GeV and pf'^'^^ 4.6 GeV. 



The possible additional bias in our used data-sets caused by the online trigger requirement 
can be assessed in the near future by analyzing the full available year 7 RHIC run MB Au+Au 
statistics recorded by STAR (see also [7 ). 

Another possibility to study the trigger bias and to select a trigger bias free jet population 
utilizing the current HT Au+Au data-set is to select the (unbiased) recoil jet. In Fig. [8] the 
di-jet A(j) distribution for Pt^f^t> 20 GeV in 0-20% central HT Au+Au and HT p+p collisions 
using the LOCone algorithm with Rc= 0.4, pf^^= 2 GeV and p^^^d_ ^ q GeV is shown. Sys- 
tematic studies are needed to utilize the di-jet pairs for further physics analysis. But given the 
current di-jet statistics it might be possible to extract a trigger bias free fragmentation function 
in Au+Au collisions. In addition di-jets in p+p collisions can be used to estimate the jet energy 
resolution in a model independent way (see for example [14J). 

In the future further systematic studies utilizing the currently available jet finding algorithms 
are needed, and will be carried out, to address systematic and conceptual uncertainties from 
full jet reconstruction in heavy-ion collisions. 
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